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Abstract 

Using 4.8 fb^ 1 of data taken with the CLEO II detector, the branching 
fraction for the Cabibbo suppressed decay D + — > tt q £ + v measured relative to 
the Cabibbo favored decay D+ -» K°l+v is found to be 0.046 ±0.014 ±0.017. 
Using V cs and V c d from unitarity constraints, we determine |/+(0)//^ : (0)| 2 = 
0.9 ± 0.3 ± 0.3 We also present a 90% confidence level upper limit for the 
branching ratio of the decay D + — > rje + v e relative to that for D + — > 7r°e + ^ e 
of 1.5. 
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Interpretation of semileptonic decays of charm mesons is theoretically straightforward. 
Amplitudes of decay modes are proportional to the CKM matrix elements and the form 
factors describing the strong interactions between the final state quarks. In this paper we 
study the Cabibbo suppressed decays D + — > 7r°£ + z/ and D + — > r]£ + v by measuring the ratios 
R w = B(D+ -> ir°l + is)/B(D + -> i?°£+z/) and i?„ = B(D + -> r]£ + u)/B{D + -> 7r°£+z/). 
Throughout this paper charge conjugate states are implied. 

The ratio i?^ is proportional to the product |/+(0)//:^(0)| 2 |14d/K s | 2 where f+(q 2 ) is the 
hadronic form factor for the decay into the pseudoscalar P. Since unitarity constraints on the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix Q yield and V^ s with good precision, mea- 
surements of R n and R v mainly provide constraints on the form factor ratios \f+(0)/ f+(0)\ 
and |/+(0)//+(0)|. Model predictions for the first form factor ratio range from 0.7 to 
1.4. This range demonstrates the difficulty in predicting how the D meson will couple to the 
light quark combinations. Assuming a monopole form |3j for the q 2 dependence of the form 
factor, the decay rate for the decay D + — > P£ + v can be written as 

V = c\\VMfm? j TT^T 2 d q \ (1) 

I 1 m*> 

where q 2 is the hadronic four momentum transfer. The mass of the nearest vector pole is 
M* = M D * for D + -> n°£u and D+ -> r/£z/, and M* = M fl , for L>+ -> 7f%/. The factor c 2 
accounts for the c?g content of the final state meson P, and is 1/2 for the tt° and r\ modes 
(dd), and 1 for the if mode (ds). There are several models that predict these rates [f|||. 
Using the framework of Heavy Quark Effective Theory and symmetry arguments, measured 
form factors from semileptonic charm decays can be compared to those for the appropriate 
b — ► u decays [fj used to extract \V u b/V c b\. 

While the Cabibbo-favored modes in charm semileptonic decay have been well measured 
lHljTj, there are relatively few measurements of Cabibbo-suppressed semileptonic decays. Pre- 
vious CLEO results for the ratio R n || are based on a total luminosity of 2.1 fb _1 , and 
are superceded by the results presented in this paper. The ratio of branching fractions 
R_ = B(D° -> tt-£ + u)/B(D° -> K~£+u) is related to R T by isospin (R w = 0.5i2_). Mark 
III 0, Fermilab E687 0, and CLEO Q have reported results for B(D° -> tt-£+z/) giving 



a current world average for i?_ = 0.102±o;[ng. 

The data sample used for this analysis was recorded with the CLEO-II detector |I2 
operating at the CESR storage ring at Cornell University. A total luminosity of 4.8 fb" 1 of 
e + e~ collisions was recorded at the T(4S') resonance and in the continuum nearby. 

In D + decays, the combinatoric background can be suppressed by requiring that the D + 
be produced in the decay chain D* + — > D + tt°. The CLEO-II detector, with its excellent 
photon detection efficiency, is ideally suited for detecting the neutral pions from this decay. 
Because the final state neutrino is not detected in semileptonic decays, we define 5m = 
M 7T °h F t+ —M hF £+, where hp refers to the D + daughter meson, the "fast" 7r° (tt f ), the K°, or 
the T). The Tig refers to the "slow" 7r° from the D* + , which is constrained by the production 
and decay kinematics to have a momentum less than 0.4 GeV/c. While the peak in 8m is not 
as narrow as the peak in fully reconstructed hadronic D + decays, a definite peak remains. 
The width of the peak in this distribution increases as more energy is carried by the neutrino. 
We therefore limit the neutrino energy by requiring 1.4 < Mh F £+ < 1.8 GeV/c 2 . 
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Electrons with momenta above 0.7 GeV/c are identified by requiring that the ratio of 
the energy (E) deposited in the Csl calorimeter and the momentum (p) measured in the 
tracking system, E/p, be close to unity and that the energy loss measured by the tracking 
system be consistent with the electron hypothesis. Muons with momenta above 1.4 GeV/c 
are identified by their ability to penetrate five nuclear interaction lengths. Electrons (muons) 
within the fiducial volume are identified with an efficiency of 94% (93%). The probability of 
a hadron being misidentified as a lepton is (0.20 ± 0.06)% for electrons and (1.4 ± 0.2)% for 
muons. We require the leptons to be found in the central region of the detector, where the 
resolution is best and the acceptance well-understood. 

Isolated photons detected by the Csl calorimeter with a minimum energy of 30 MeV 
are paired to form ir° and rj candidates. For the slow pion, the 77 mass is constrained to 
be within 2.5 standard deviations (about 12.5 MeV/c 2 ) of the nominal 7r° mass. For the 
fast 7T° (77), the reconstructed mass is required to be within the range 0.105-0.165 GeV/c 2 
(0.510-0.585 GeV/c 2 ). The decay channel 77 — > 7r + 7r~7r° was not considered because of its 
low reconstruction efficiency. For the normalizing D + — > K^i^v mode, we identify the K° 
through the 7r + 7r~ decay of its K$ component. We require the 7r + 7r _ pair to form a secondary 
vertex of the correct mass that is displaced at least four standard deviations from the primary 
vertex. 

Combinatoric backgrounds are reduced by several means. We impose the kinematic crite- 
ria 0.175 < p n o < 0.350 GeV/c, p hp > 0.7 GeV/c, and \p hp + p e \ > 2.1 GeV/c. Backgrounds 
from B meson decay are reduced by requiring that the ratio of Fox- Wolfram moments [ 13 
i?2 = H2/H0 satisfy R2 > 0.2. Finally, we consider only well-measured tracks and events 
with a hadronic event structure. 

Backgrounds can be divided into four classes: fake slow pions (fake D*s), fake fast 
hadrons, fake leptons, and uncorrelated fast-hadron, lepton pairs (fake D + s). The major 
contribution to the fake D + background in the D + — > 7r°£ + z/ channel comes from feed- 
down from D + — > K°l + v, K° — > 7r°7r°. We can correct for this background knowing only 
the ratio of the reconstruction efficiency for D + — > tt°£ + u to the efficiency to reconstruct 
D + — > K°£ + u, K° — > 7T 7r° as 7r°£ + z/, which we determine from Monte Carlo simulation. 
Monte Carlo studies indicated that the feedthrough from other semileptonic charm decays 
and from BB events is negligible. The other background components were determined from 
fits to the data. 

We fitted the two dimensional distribution of 5m versus fast hadron mass to extract 
the signal yield. Figure |I| shows the distributions for the signal Monte Carlo and data for 
the D + — > 7T°e + z/ e mode. Figure § shows the 5m projection for the ir°, K°, and 77 modes. 
Figure shows the fast hadron mass distributions. The fits used a parametrization of the 
fast hadron mass obtained by fitting these one- dimensional projections. The signal shape in 
5m was determined from fits to the distributions of reconstructed signal Monte Carlo. The 
fake lepton background was determined by performing a fit to the distributions of events 
which satisfied all requirements except for the lepton identification requirement. The signal 
yields from these fits were then scaled by the measured misidentification probabilities and 
subtracted from the yields from the fit to the data. The parameterization of the fake D* 
background in 5m was determined by looking at a sample of data events whose fast hadron 
mass was more than 4 sigma from the nominal mass. The signal yields, fake lepton yields, and 
signal reconstruction efficiencies are presented in Table |. The efficiencies were determined 



5 



from fits to the distributions from samples of reconstructed signal Monte Carlo. 

With the results from the fits given in Table |, we proceed to calculate the ratio of 
branching fractions R n = [B(D + — > 7t°£ + i>)}/{B(D + — > K°£ + v)]. For each leptonic mode we 
define, 

N(ti f £ + u) e^ e+u (W£+v) e^ + >°l+z/) 
* N(K° s £+u) e{7r°£+u) e(n°£+v) 

Here N{tx f £ + v) and N(Kg£ + u) are the two signal yields after background subtraction, 
e{ix Q £ + v) is the efficiency for a 7r°£ + i/ decay to be reconstructed as itself, ej^o £+v (K°£ + u) 
is the efficiency for a K £ + u decay to be reconstructed as K°£ + i>, and e-^ +[y (7r°£ + z/) is the 
efficiency for a K°£ + u decay to be reconstructed as n°£ + v. The ratio for electrons was found 
to be R n = (4.5 ±1.6 ±1.9)%, where the first error is statistical and the second is systematic. 
The ratio for muons was found to be = (4.8 ±3.1 ±3.2)%. Here the error from fake muon 
subtraction is substantial and the detection efficiency is lower than for the electron channel. 
We combine the results weighted by their errors to find R n = (4.6 ± 1.4 ± 1.7)%. 

Most of the systematic effects cancel in the ratio of branching fractions because we 
impose similar requirements on both the signal and normalization modes. The systematic 
error for the electron channel is dominated by the parameterizations of the shapes in the 8m 
distribution (30%). This error is correlated between the ir°e + v e and K°e + v e channels. The 
systematic error in the ratio due to Monte Carlo simulations of K° s — > 7r + 7r~ and tt f — > 77 
is conservatively placed at 10%. Other systematic errors for the electron channel include: 
statistical error on efficiency fits from Monte Carlo samples (7%), fake lepton subtraction 
(7%), D + — > K°e + u e feeddown (9%), other semileptonic charm decay feeddown (16%), and 
BB feeddown (13%). The systematic errors are added in quadrature to obtain a total 
systematic error in the ratio for electrons of 41%. 

The fit to the D + — > r]e + v channel yielded 6 ± 8 events. We did not consider the muon 
channel due to the low detection efficiency. To obtain an upper limit on R v , we scale this yield 
by the reconstruction efficiency of (0.26 ±0.02)%, and normalize to the average D + — > tt°£ + v 
yield of (4.39 ± 2.22) x 10 3 events. The latter was estimated from our R v measurement and 
the average of the efficiency-corrected yields for D + — > K°£ + i> in the electron and muon 
channels. We find R v = B(D+^nO£+u) =< ^ a ^ ^he ^0% confidence level. This result is 
dominated by statistical error, but includes a 30% systematic error that was combined in 
quadrature with the statistical error. 

We have measured the branching fraction of the Cabibbo suppressed decay D + — > 7i £ + u 
relative to D + — > K°£ + u. Using our measurement of this ratio, we find using Equation 
(1) |^(0)//^(0)| 2 |V; d /1/ cs | 2 = 0.046 ± 0.014 ± 0.017. The integral in Equation (1) times 
the constant term is approximately 1 here. Unitarity constraints on the CKM matrix yield 
\V cd /V cs \ 2 = 0.051 ±0.001 0] which translates to a value of 0.9 ±0.3 ±0.3 for |/^(0)//^(0)| 2 . 
Model predictions || are in agreement with our measurement. We can combine our mea- 
surement of R n with the measurements of 0.5 x R_ to obtain R n = 0.050 ± 0.008 and 
= 0-99 i 0.08. The upper limit on the ratio R^ is consistent with current 

predictions. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. J.P.A., J.R.P., and I.P.J.S. thank the NYI program of 
the NSF, M.S. thanks the PFF program of the NSF, G.E. thanks the Heisenberg Foundation, 
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TABLE I. Results of fits to the M^ F versus 5m distributions for each of the three analyses. 



Sample 












Electrons 


Muons 


Electrons 


Muons 


Data 


75 ± 15 


83 ±18 


530 ± 29 


178 ± 17 


n a kp T iPntnn 

A ClIVC UCUIU11 




48 + 1 


7+2 


25 + 5 


SUBTRACTED 


65 ± 15 ± 20 


35 ± 18 ± 16 


523 ± 29 ± 38 


153 ±17 ±13 


e{7r°e+i> MC)% 


1.01 ±0.05 ±0.03 


0.66 ±0.05 ±0.01 






e{K®l + v MC)% 


0.020 ± 0.004 


0.007 ±0.005 


0.54 ±0.01 ±0.02 


0.17 ±0.02 ±0.01 


e{K*°£ + v MC)% 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


YIELD (xlO 3 ) 


6.44 ± 1.49 ±2.00 


5.30 ± 2.73 ± 2.84 


96.85 ± 5.37 ± 7.95 


90.00 ± 10.00 ± 13.66 
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and V.S. thank the A. P. Sloan Foundation, and A.W. and R.W. thank the Alexander von 
Humboldt Stiftung for support. M.S. is supported as a Cottrell Scholar of Research Corpo- 
ration. This work was supported by the National Science Foundation, the U.S. Department 
of Energy, and the Natural Sciences and Engineering Research Council of Canada. 
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FIG. 1. The distribution of M 77 versus 5m for a) D + — > ir°e + is Monte Carlo events and b) data. 



REFERENCES 

[1] R.M. Barnett et al. (Particle Data Group), Phys. Rev. D54, 1 (1996). 

[2] S. Narison et al, Phys. Lett. B377, 163 (1994); V. Lubicz et al, Phys. Lett. B274, 415 
(1992); M. Bauer, B. Stech, M. Wirbel, Z. Phys. C34, 103 (1987); C.A. Domingues, 
Phys. Lett. B207, 499 (1988); M. Crisafulli et al, Phys. Lett. B223, 90 (1989); G.P 
Lepage and S.J. Brodsky, Phys. Rev. D22, 2157 (1980). 

[3] A. Ali and T.C. Yang, Phys. Lett. B65, 275 (1976). 

[4] J.G. KSrner, K. Schilcher, M. Wirbel, and Y.L. Wu, Z. Phys. C38, 511 (1988). 
[5] N. Isgur and D. Scora, Phys. Rev. D52, 2783 (1995); see also N. Isgur et al, Phys. Rev. 
D39, 799 (1989). 

[6] N. Isgur and M.B. Wise, Phys. Lett. B232, 113 (1989); N. Isgur and M.B. Wise, Phys. 
Rev. D42, 2388 (1990). 

[7] A. Bean et al, Phys. Lett. B317, 647 (1993). 

[8] M.S. Alam et al, Phys. Rev. Lett. 71, 1311 (1993). 

[9] J. Adler et al, Phys. Rev. Lett. 62, 1821 (1989). 
[10] PL. Frabetti et al, Phys. Lett. B382, 312 (1996). 
[11] F. Butler et al, Phys. Rev. D52, 2656 (1995). 
[12] Y. Kubota et al, Nucl. Intr. and Meth. A320, 66 (1992). 
[13] G.C. Fox and S. Wolfram, Phys. Rev. Lett. 41, 1581 (1978). 



8 



3030197-002 




0.14 



0.26 



0.30 



0.18 0.22 

<5M (GeV / c 2 ) 

FIG. 2. The 5m spectra for data passed through the a) D + 
0.115 < M n o < 0.153 GeV/c 2 , b) D+ -► We + u e analysis with 0.48 < M K o 



7r°e- 



> e analysis with 
< 0.52 GeV/c 2 , 



and c) D + -> r/e+^ e analysis with 0.51 < M v < 0.58 GeV/c 2 . 
while the dashed line indicates the background function. 



The solid line indicates the total fit 
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FIG. 3. The a) M 77 = M^o , b) M^+^-, and c) M 77 = M rj spectra for data passed through the 

D + — > 7r°e + u e , D + — > K°e + f e , and Z) + — > r]e + u e analyses, respectively, with (5m < 0.3 GeV/c 2 . 
The solid line indicates the total fit while the dashed line indicates the background function. 
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